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Abstract: We study the wave equation of a massive vector boson in the background of 
a L)-dimensional Schwarzschild black hole. The mass term introduces a coupling between 
two physical degrees of freedom of the field, and we solve the resulting system of ODEs 
numerically, without decoupling. We show how to define decoupled transmission factors 
from an S-matrix and compute them for various modes, masses and space-time dimensions. 
The mass term lifts the degeneracy between transverse modes, in D = 4, and excites the 
longitudinal modes, in particular the s-wave. Moreover, it increases the contribution of 
waves with larger £, which can be dominant at intermediate energies. The transmission 
factors are then used to obtain the Hawking fluxes in this channel. Our results alert for 
the importance of modelling the longitudinal modes correctly, instead of treating them 
as decoupled scalars as in current black hole event generators; thus they can be used to 
improve such generators for phenomenological studies of TeV gravity scenarios. 
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1 Introduction 

Separation of variables and decoupling of degrees of freedom are two key properties in the 
study of wave equations in gravitational backgrounds. Most remarkably, the gravitational 
perturbations of the Kerr black hole both separate and decouple, as first exhibited in 
the celebrated work of Teukolsky [1], which allowed to show linear stability of the Kerr 
solution. Introducing charge or considering higher dimensions, for instance, seems to spoil 
such exceptionality; indeed gravitational perturbations do not seem separable in either 
the four dimensional Kerr-Newman background [2] or the higher dimensional Myers-Perry 
background [3]. 

Failure to achieve separation of variables or to decouple individual degrees of freedom 
typically leads to the problem of solving sets of coupled ordinary differential equations 
(ODEs).^ In the absence of analytic strategies, a full solution of the problem can only 
be obtained numerically. Surprisingly, numerical solutions for such (common) problems 
seem to be essentially unexplored, in the context of Hawking radiation, in contrast to the 
vast literature concerning fully decoupled and separable problems. Perhaps the reason is 
connected to the problem of quantisation, which is usually straightforward when a complete 
set of decoupled modes is available. We will show, however, that coupled systems (of the 
type considered herein) may be treated similarly using an S-matrix type formalism which 
allows decoupling in the asymptotic regions. 

""^ Non-separability leads to partial differential equations but these may be discretised to systems of ODEs. 
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In this paper we shall consider a wave equation which separates but does not decouple, 
in the background of a D-dimensional spherically symmetric Schwarzschild black hole [4]: 
the massive spin 1, or Proca, field. In contrast, observe that the Maxwell field in the 
D = 4 Schwarzschild background has two physical decoupled polarisations. Introducing 
the mass term, the spin 1 field gains a longitudinal polarisation, and two of the physical 
degrees of freedom cannot be decoupled. In higher dimensions the situation is similar, 
when set in appropriate variables. We shall solve the wave equation for the Proca field in 
the Z)-dimensional Schwarzschild background numerically to obtain a scattering matrix. 
This will allow us to define transmission factors as well as the Hawking radiation fiux of 
vector bosons from such black holes. 

The technique used herein can in principle be applied to other problems of test field 
wave equations or studies of gravitational perturbations, including the computation of 
quasi- normal modes^. Our choice of the Proca field was primarily due to its simplicity, but 
most importantly for its phenomenological interest to TeV scale gravity scenarios [7-10] . In 
such scenarios, scattering processes with centre of mass energy well above the fundamental 
Planck scale, should be dominated by classical gravitational interactions [11]. Moreover, 
for sufficiently small impact parameter black holes should form in particle collisions, as 
made manifest by numerical evidence [12]. This was predicted by early trapped surface 
calculations [13] in D > 4, which gives bounds on the gravitational radiation emitted; the 
latter was also estimated through perturbation theory studies of head on collisions of shock 
waves in = 4. In D > 4, recent results have produced improved semi-analytic estimates 
within perturbation theory [14] . Other estimates may be obtained from high energy black 
hole collisions in numerical relativity. The first results for numerical collisions of higher 
dimensional black holes have been released in [15-17], based on the frameworks in [18, 19]. 
These black holes would then decay via Hawking radiation [20] and leave observable signa- 
tures [21], such as large multiplicity of jets or large transverse momentum, and considerable 
missing energy due to gravitational wave emission. To better understand and model these 
scenarios, as well as to improve the phenomenology of black hole event generators such as 
CHARYBDIS2 [22] and blackmax [23], currently in use at the Large Hadron Colhder, it is 
important to understand the Hawking radiation from higher dimensional black holes. This 
is especially relevant because the first experimental tests of trans-Planckian gravity sce- 
narios are under way [24-26] , and therefore improving their phenomenology is very timely. 
The computation of Hawking radiation for such models requires solving the wave equation 
for various fields with different spins, masses and charges in these black hole backgrounds. 
In particular, the Proca field in the D-dimensional Schwarzschild background has not been 
considered up to now. 

Our results exhibit distinctive features as we introduce the mass term, such as the 
lifting of the degeneracy of the two transverse modes in four dimensions, the appearance of 
longitudinal mode contributions (absent for Maxwell's theory) and in particular the s-wave. 
One feature that appears not to have been discussed in the literature is that in four and 
five space-time dimensions, the transmission factor has a non- vanishing value in the limit 

^See [5, 6] for studies of quasi-normal modes where similar coupled systems have been considered. 
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of zero spatial momentum. We will also find the expected suppression with M, but perhaps 
the most relevant feature is to notice the increasing importance of the longitudinal modes 
and larger i partial waves. Moreover, the precise contribution of the longitudinal modes 
deviates from simple models that have been used in the aforementioned Monte Carlo event 
generators to produce phenomenological data, providing our results with phenomenological 
interest in this context. 

This paper is organised as follows. In Section 2 we explain how to decompose the 
Proca field equation in a generic factorizable background geometry consisting of a warped 
product of an m-dimensional space with an n-dimensional Einstein space, using the har- 
monic functions of the Einstein space. In Section 3 we specialise the equations of motion to 
a D-dimensional Schwarzschild background and present an appropriate set of independent 
degrees of freedom of the Proca field. We then study the asymptotic and near horizon 
behaviour of the coupled radial equations, which is required to impose the boundary con- 
ditions and to define the scattering matrix. In Section 4 we discuss how the scattering 
matrix is used to compute the transmission factor and the Hawking spectrum. In Section 5 
we discuss the numerical method and results, and we conclude with a discussion in Section 
6. Some technical relations are left to the Appendix. 

2 The Proca equations in Einstein symmetric spacetimes 

In this Section we present the wave equations for a Proca field which may be complex and 
charged under a U{1) electromagnetic field. This covers the effective fields describing the 
Z and W particles in the Standard Model of particle physics (SM), the former being a 
neutral Proca field and the latter being an electromagnetically charged Proca field. In the 
next Section we will specialise to the neutral case which we want to study in detail. The 
Lagrangian is 

C = -^^Wl^W^" + M'^WpV^' + iQWlW^F^''' , (2.1) 

where Wf^i, = d^W^ — d^W^ and we have included the coupling of the W to the electro- 
magnetic field strength tensor F'"^ as in the SM^ . The equations of motion for W when all 
the background fields are fixed are 

V^Wf"" - M^W^" - iQW^F^"" = . (2.2) 

For the gravitational background, in this Section, we consider Einstein symmetric spaces 
of the form [27] 

ds^ = hMdy'^dy' - r{yfdal , (2.3) 
where cj„ is an n-dimensional Einstein space with constant curvature K and 

dal = aij{x)dx^dx^ . (2.4) 

Later we will specialise to the Schwarzschild black hole. We use indices a,b,c, . . . for 
the first set of coordinates, {y""}, spanning the m-dimensional space with metric hab{y)] 

^In this paper we use the particle physics convention for the signature of the metric (H . . .). 
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and indices k, . . . for the second set of coordinates, {x*}, spanning the Einstein space. 
Furthermore, geometric quantities and differential operators on an are denoted with hats. 
This covers several interesting cases such as 2 + n-dimensional spherically symmetric black 
holes or a singly rotating black hole in 4 + n-dimensions (Kerr, Myers- Perry, etc . . . ). 

To write down the equation of motion we use a decomposition of the vector field in 
tensorial types [27]. Wa are m-scalars, with respect to (7„, so they must obey 

(A + 4)Wa = 0, (2.5) 

(kq is the spin-0 eigenvalue). Wi is a co- vector field which can be decomposed into a scalar 
<I>, and a transverse co- vector , i.e. 

Wi = Di<P + , DiW'^' = , (2.6) 

where Di is the covariant derivative on cr„ and we use an explicit hat to denote raising of 
indices using aij, i.e. note that = —r{y)'^W'^^. Since $ is a scalar, it obeys (2.5). The 
transverse vector obeys 

[a + kI)w[ = 0, (2.7) 

where now is the spin-1 eigenvalue. This decomposition allows for an expansion of 
the various degrees of freedom Wa,^,Wf in a basis of harmonics of the Einstein space. 
Furthermore, this decomposition allows for a decoupling of the field equations into an 
independent vector mode Wif and m + 1 coupled scalar fields for each set of quantum 
numbers labeling the basis of harmonic functions. Observe, however, that not all these 
modes correspond to physically independent degrees of freedom, as shown below. 



3 Neutral Proca field in a spherically symmetric black hole 

Let us now specialise to the Q = case. We shall consider separately two cases according 
to the value of kq. 

3.1 Modes with kq ^ 

Expanding the field equations with the decomposition (2.6), using conditions (2.5) and (2.7), 
and defining Ba = Wa — we obtain 



^0 R ^af o 
-^Ba —Ob 



+ M^Ba + da 



db ( r^-^h'^B, 



, (3.1) 



1 



R 



r.n—2 



1 



da ( r^'-^^^Bt, ) - M''^ = , (3.2) 



n 



W^^ = . (3.3) 



We consider the spherically symmetric case with {y"} = {t,r}, \h\ = 1, hat is diagonal, 
htt = —1/hrr = V = 1 — n/r'^~^ and ^ = r^~^. We choose units such that the horizon 
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radius is rn = 1- Since <I> is given by the second equation in terms of the other fields, it 
is a non-dynamical degree of freedom. In four dimensions, this agrees with the fact that 
a spin-1 massive field has three possible physical polarizations which in this case will be 
the two dynamical scalars and the transverse vector. In higher dimensions, the transverse 
vector on the n-sphere will contain more (degenerate) polarizations. 
We can factor out the spherical harmonics through the expansion 



W 



Ti 



/3^(y)3^A(x) 



(3.4) 



where A denotes the mode eigenvalues for the corresponding harmonic functions. Further- 
more, making the ansatz 



A 



-iojt 



and using (3.1) and (3.3) we find 

1 d 

dr 



dr 



n-2 



-n—2 



-iut 



V 



V 



r"- dr \ dr , 



-^ + M'\V 



r 



= e-^'"*T(r) 
X — itoV'ip 



V d 



>.n— 2 



dr 



,n-2 



V- 



dr 



1 n 



+ V 



T = 



(3.5) 
(3.6) 

(3.7) 



Thus we obtain two second order coupled radial equations for {-0, %} and a decoupled 
equation for T. Note that Kq = £(£ + n — 1) and n\ = l{t + n — 1) — 1 with £ starting at 
zero and one respectively. The third combination is + = i{l + n — 1) + n — 2. 

The manipulations leading to the two coupled equations above are only valid for non- 
zero M. In the exactly massless (Maxwell) theory, a similar calculation leads to a single 
decoupled equation for one of the scalar modes which is 



V d 

.n-2 ^j. 



^n-2y 



(3.8) 



whereas the other mode V = iVd^ 



^x)/('^'"" ^) is non-dynamical. Here dr = d/dr. 



The transverse mode - ruled by equation (3.7) - remains the same for any M; in particular, 
for M = 0, and (only) n = 2 it becomes equivalent to (3.8). This will be manifest in the 
numerical results. 



3.2 Modes with kq = 

For the exceptional modes with kq = 0, $ does not enter the wave equation so it is a free 
non-dynamical field. The corresponding equation for wj^^ is (the superscript denotes it is 
the exceptional mode) 



haf 



(3.9) 
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When 7^ one uses an ansatz similar to the previous section to obtain a radial equation 
for a dynamical degree of freedom 



V d 
r" dr 



d 



2 - M2y dr 



+ 1 



(0) 



0, 



(3.10) 



and a non-dynamical one, x''^'* = iu^Vjiuj^ — M^l/)^,.^^*'). Otherwise, for = 0, we 
recover the well known result that all the exceptional modes are non-dynamical (see e.g. 
[28]). 

Now that we have covered all possibilities, several comments are in order. Firstly there 
is a discrete difference between the small mass limit and the exactly massless theory since we 
have different sets of equations for each case. This should not be surprising since there is an 
extra longitudinal mode for massive vector bosons. Secondly, the equations for the Maxwell 
theory case are all decoupled, in agreement with previous work [29]. Since decoupled radial 
equations have been extensively studied in the literature we will not present the details of 
our analysis of such modes and refer to the method in [30, 31]. Thus in the remainder we 
will focus on the solution of the coupled system for the massive theory, which will be used 
in conjunction with the decoupled modes to obtain the full Hawking spectrum in Section 
5. 

3.3 Boundary conditions and radial system 

In this section, we start by finding a series expansion of the solution near the horizon for 
the coupled system {V^,x}- This will be used to initialise the corresponding fields for the 
radial integration at r = 1.001. If we define ?/ = r — 1, Eqs. (3.5) and (3.6) become 



d^ d 
M(T)—^^N{r)—^P(T) 
dy dy 



d^ 



d 



M{r)—. + N{r)—+P{r) 
dy dy 



^p + Q{r)x = , 
X + g(r)V' = , 



(3.11) 
(3.12) 



where the polynomials are defined in the Appendix. Making use of Frobenius' method to 
expand -0 and x ^ 

00 00 

4^ = y'^^f^jy^ , x = y''^^jy^ , (3-13) 

j=0 j=0 

and inserting the above two equations into Eqs. (3.11) and (3.12), we obtain 



lUJ 



or p=l± . (3.14) 

1 n — 1 

We want to impose an ingoing boundary condition at the horizon, so we must choose the 
minus sign. Furthermore, after this sign choice, the right hand side case produces a series 
expansion which is a special case of the left hand side (where the first coefficient is set to 
zero), so without loss of generality we choose p = —ioj/{n — 1). One then writes down 
the recurrence relations and conclude that a general solution close to the horizon can be 
parametrised by the two coefficients and i^i. The other coefficients are generated by the 
recurrence relations (A.l). 



P 



n 
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To understand the asymptotic behaviour of the waves at infinity we now study a large 
r asymptotic expansion in the form 



j=0 



X 



j=0 



Inserting this into Eqs. (3.5) and (3.6) we obtain, at leading order, 



/3 = ±ik 



1 



n 



ziziip or p 



n 



(3.15) 



(3.16) 



where ip = 5n.2{^'^ + ^^)/(2A;). Thus one can show that asymptotically^ 



r 2 



an + 



+ 



+ 



r 



X 



r 2 



k c+ 
a; r 



+ . . . 



A; c 



+ ... 



(3.17) 



where $ = /cr + (/?logr and is defined in the Appendix, Eq. (A. 2). So as expected, 
each field is a combination of ingoing and outgoing waves at infinity. This asymptotic 
expansion also shows that for a generic wave at infinity, we can choose four independent 
quantities {a^, a^} at infinity, to characterise the solution. This is expected, since we have 
two coupled scalar fields and for each scalar degree of freedom we must have an associated 
ingoing wave and outgoing wave. Thus we can define four new fields {x^jV'^} (which 
will asymptote respectively to {a^,a]^}), by truncating the expansion for the fields and 
for their first derivatives at infinity. Such a transformation can be written in matrix form 
by defining the 4- vector = {ipj^,'ip-,x+^X-) for the new fields, and another 4- vector 
V"^ = (^, dyip, X) drX) for the original fields and derivatives. Then the transformation is 
given in terms of an r-dependent matrix T defined through 



V = T* 



(3.18) 



which we provide in the Appendix, Eq. (A. 3). Finally we obtain a first order system of 
ODEs for the new fields. First we define a matrix X through 

f=XV, (3.19) 

which is read out from the original system (3.5), (3.6). Its explicit form is given in the 
Appendix, Eq. (A. 4). Then we obtain 



dr 



dr I 



(3.20) 



We can write other equivalent systems using different T matrices. In particular, we have 
also integrated a first order system using the fields ijjs = kip — isdrtp and Xs = kx — isd^x 
which produced numerically equivalent results. The only difference is that for such fields 
we need to extract 0(r~^) coefficients to obtain af . 



*We have used, without loss of generality, the leading power behaviour for p and discarded the second 
option similarly to the near horizon expansion. 
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4 The Hawking spectrum 



The boundary conditions we have chosen in Section 3.3, are suitable for the computation of 
the Hawking spectrum of radiated quanta from the black hole. The Hawking spectrum is 
generically given by a sum over a complete set of modes with labels of the transmission 
factor times a thermal average number of quanta produced at the horizon (n^). This is 
defined for a basis of decoupled modes. In our problem, we have a sub-set of modes, the 
transverse vector mode, and the i = {kq = 0) mode, which are decoupled. But we also 
have a tower of modes which are coupled two by two for each i > 0, the two scalars ip and 
X- It is not obvious how to decouple them for all r through an explicit transformation. 
Instead, let us try to understand how to extract the relevant information in the asymptotic 
regions. 

Let us denote the two coupled fields by a 2-vector U'^ = (ip, x) and represent the cou- 
pled system of radial equations through a (linear) second order matrix differential operator 
acting on U, i.e. D^^^U = 0. The system is coupled because of the off diagonal ele- 
ments of the P^2^ operator. To decouple the system we would have to find a transformation 
of the fields U = ^U, such that the new differential operator P^^^ = D^'^^oA is diagonal. 



Even without finding such a transformation explicitly, one can draw some conclusions by 
assuming its existence^. In particular we may establish a map between our general solution 
of the coupled system and the actual decoupled solution, for each of the asymptotic regions 
(horizon and far field). To find such a map let us first summarise the information we have 
on the general solution of the coupled system. 

In Section 3 we have found that a general solution is parametrised by 4 independent 
coefficients in one of the asymptotic regions; either at the horizon or at infinity. Once we 
have chosen one set of coefficients, say at the horizon, due to the linearity of the equations, 
the 4 independent wave components at infinity are a linear combination of the 4 coefficients 
at the horizon. Let us formally denote the ingoing and outgoing wave coefficients at the 
horizon (+/— respectively) by 



where i = 1,2 since we have two fields. Similarly, the coefficients at infinity are defined 
as the large r limit of the ^ field components (up to linear transformation which we will 
define next), i.e. 



^In fact, for example, if we consider ^ to be a general r-dependent matrix, we can write down two 
conditions for the four arbitrary functions of such a matrix. Thus, in principle, there is enough freedom. 



I.e. 




(4.1) 



h = (h+,h-) = (/i+,/i-) 



y = (y^,y ) = {yl,yi ) 
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with i = 1,2 for ip and x respectively. Due to linearity, we can define a scattering matrix 



Sh ^ 






-+ 


s+- 




+ 







(4.2) 

which is a set of numbers (depending on energy, angular momentum, etc. . . ) containing 
all the information on the scattering process. It can be fully determined by considering 
specific modes at the horizon and integrating them outwards. In our problem, we have 
imposed an ingoing boundary condition at the horizon which is simply = 0. Then 

y** = S'-h- . (4.3) 

Taking the s = — component, and denoting the inverse matrix of S by (S )~^, we 
invert (4.3) to obtain the wave at the horizon given the ingoing wave at infinity 

h~ = {S-)-'y- . (4.4) 

Inserting this relation back in the s = + component of (4.3), we obtain the outgoing wave 
in terms of the ingoing wave, at infinity 

y+ = S+-(S— )-V"^R-y" , (4.5) 

where in the last line we have defined the reflection matrix R. Before proceeding, we note 
that there is still some freedom in the definition of the asymptotic coefficients since any 
(non-singular) linear combination is equally good from the point of view of satisfying the 
boundary condition. This freedom can be written in terms of 3 matrices M**, M]^ relating 
some new fields (hatted) to the old fields 

y^ = M^r , h- = M^h- . (4.6) 

Since this represents the most general parametrisation of the solution in the asymptotic 
regions, there must be a choice which decouples the fields in those regions. To find the 
correct transformation we need a physical prescription. 

To obtain the transmission factor for the decoupled components, it is instructive to 
remind ourselves of the calculation of the transmission factor for a single decoupled field. 
It is defined as the fraction of the incident wave which is transmitted to the horizon. If we 
look at a wave with energy uj (for an observer at infinity), with ingoing/outgoing amplitudes 
Y^°°\ then [32] 



T = ' - 11^+ I = ^ V ^ = gL , (4.7) 

where in the last step we note that T can be re-expressed as a ratio between the total 
incident energy fiux (which is the difference between the energy carried by the ingoing 
wave and the energy of the outgoing wave) and the incident energy fiux associated with 
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the ingoing wave at infinity {J-'^)- The former is the flux of energy transmitted down to 
the horizon. 

We now compute the energy fluxes through a sphere at radius r using the energy 
momentum tensor. This will allow us to identify the decoupled fields at infinity and at 
the horizon, and in particular, the ingoing and outgoing decoupled waves at infinity. Such 
a flux is shown to be conserved in our background, by using the conservation law for the 
energy momentum tensor, combined with the fact that the spatial integral of T^* for each 
energy eigen-mode is constant. It is defined, evaluated at r, as 



J'\r = - d^T/ (4.8) 

where is the volume element on a t, r = constant hyper-surface; up to an irrelevant 
normalisation, the energy momentum tensor for the complex neutral Proca field is 

r^^ = -- (^lytM"!^^ _ mV'W + c.c.j - . (4.9) 



If we insert this in (4.8), assume a field configuration with a well defined energy uj, and 

.2 
^0' 



make use of the equations of motion, then, for the non-trivial case of / 7^ Kq, we 



obtain 

A A ^ 



2 



+ c.c 



2M2 ^ 

(4.10) 

where for convenience we define ^ = dril^ + iwx- Modes with different angular momentum 
eigenvalues are clearly decoupled, as are the transverse vector mode contributions in the 
first sum. The terms in the second sum couple two fields for fixed A. We can compute the 
fiux at infinity and close to the horizon and express it in terms of the asymptotic coefficients 
in the corresponding region. Focusing on a specific mode and in the coupled part of the 
fiux (second sum in (4.10)) 

-^r^'^" = I2/0 1' - \yt? + \yi? - \yt? = (y")V- - (y+)V+ , (4.ii) 

where yf are linear combinations of the asymptotic coefficients a| given in the Appendix, 
Eq. (A. 5). This choice of y| is already in a form close to decoupled, since we have separated 
the modulus square of the incident contribution from the reflected contribution, without 
interference terms. This form is invariant under separate unitary transformation of y^. 
Using the reflection matrix we obtain 

^coupled ^ (y-)t (^1 _ RtR^ y- ^ (y-)tTy- , (4.12) 

where we have defined a (hermitian) transmission matrix T. This can be diagonalised 
through a unitary transformation which is the remaining freedom we have for y~. In fact we 
can do even better, and diagonalise the reflection matrix R with a bi-unitary transformation 
using the arbitrary unitary transformations. Then the fields are manifestly decoupled 
at infinity, both at the level of the reflection matrix and the transmission matrix. As 
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a consequence, in the decoupled basis, an incident wave is reflected back in the same 
decoupled mode without interference with the other mode. Finally, the transmission factors 
are simply the eigenvalues of T, since they are each associated with a decoupled component. 

Furthermore, one can use the conservation law for the flux, to find an alternative ex- 
pression for the transmission matrix, at the horizon (this will be useful to control numerical 
errors). The total flux at the horizon is 

J^°"P'"^ = (h-)^h- , (4.13) 

where the coefficients are linear combinations of the two independent Ui coefficients 
{i = 0,1), given in the Appendix, Eqs. (A. 6). Eq. (4.13) establishes the important point 
that the flux is positive definite, so the transmission factors must be positive definite (as 
expected since there is no superradiance in Schwarzschild spacetime). Finally, using the 
relation between y~ and h~ through^ S~~, we find 

T = (S—S^—)-^ . (4.14) 

Once we have obtained the transmission factors, the number and energy fluxes are given 
by the standard result 

d{N,E} 1 v-v^ {l,uj} 



^EE .J';:L/ cTc. (4.15) 



dtduj 2tt ^ ^ exp{ijj /Th) — 1 
^ C 

where C is a label running over the final set of decoupled scalar modes and the transverse 
mode, and are the degeneracies of the corresponding spherical harmonics. Labeling the 
scalar and vector harmonic degeneracies by ds and dy respectively we have [33] 

(n + 2^-l)(n + ^-2)! 

{n + 2i-l){n + e-l){n + i-3y. 

= — {e+m-mn- 2)1 — • (^-^^^ 

The Hawking temperature in horizon radius units is 

Th = ^. (4.18) 

5 Results 



In this Section we present a selection of numerical results to illustrate the behaviour of the 
transmission factors and the corresponding Hawking fluxes. To integrate the coupled and 
decoupled radial equations, we wrote firstly test codes in mathematicaT and then a code 
in the CH — h language, using the numerical integration routines of the Gnu Standard Li- 
brary (GSL). Besides using different programming frameworks we have also tested different 
integration strategies which all agreed within relative numerical errors smaller than 0.1 %. 

®Note that the relation between h~ and (M^)^^y^ can be made diagonal using M^, so the problem is 
also decoupled at the horizon. 
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In fact, most of our numerical points have a precision which is one order of magnitude 
better. To check numerical errors we have integrated the radial equations up to a large 
radius of typically r = IQ'^rn and varied this up to a factor of 3 to check the precision. 
Furthermore we have used the two expressions for the transmission factor from Eq. (4.12) 
and (4.14) which agree within the quoted precision for almost all energies. The exception 
is for small energy, where the first definition converges poorly. This can be explained by 
a simple analysis of propagation of errors combined with the fact that the coefficients 
grow very fast as we decrease energy, thus requiring a very large precision for some fine 
cancellations to occur. The second expression is thus more natural in that limit since it 
does not need such cancellations and does not require such large precision. 

We have generated several samples of transmission factors, some of which are displayed 
in Fig. 1. Hereafter, we shall denominate the partial waves associated to the different modes 
of the Proca field by ii,i2,^T and i = 0, where £±,£2 correspond to the two coupled modes 
ruled by Eqs. (3.5) and (3.6), £t to the decoupled mode described by Eq. (3.7) and i = Oto 
the kq = mode, ruled by Eq. (3.10). Moreover, partial waves associated to the Maxwell 
field shall be denoted by ^e, and are ruled by (3.8). 

In the top row panels of Fig. 1, we show the partial wave contributions for n = 2, 3, 4 
in the zero mass limit. Some general properties are as follows. The curve becomes 
shifted towards higher frequencies both as i is increased, for n fixed, and as n is increased, 
for i fixed. The former can be understood from standard geometrical optics arguments. 
Moreover, for this choice, there is always a numerical coincidence between one of the 
partial waves (ii) obtained from the two coupled fields and the electromagnetic partial 
wave Ie- The i = and £2 modes are always absent in the Maxwell theory, so they can 
be associated with the longitudinal polarization of the massive vector field. Similarly, the 
It and ii partial waves are associated with the transverse polarizations of the field. A 
qualitative dependence on dimension is that for n = 2, £t and £1 (or £e) modes are all 
equal. Curiously, this is in agreement with the fact that they describe the same number of 
transverse degrees of freedom as can be seen from the degeneracies (4.16) specialised for 
n = 2. This degeneracy is lifted for n > 2. 

For non-zero mass (middle and bottom row panels of Fig. 1), the degeneracy observed 
for n = 2 in the massless limit is lifted. Also, we observe, for all n, that modes with 
higher £ partial waves (especially £1 modes) become a more dominant contribution at 
lower energies, as compared to lower £ partial waves of other modes. In particular for 
M = 1, the transmission factor for £1 becomes the largest for small energy. This effect of 
excitation of sub-dominant partial waves is well known to exist for example as we increase 
n (and we can also observe such effect in our plots) as well as with the introduction of 
black hole rotation [34]. If this effect persists cumulatively on a rotating background, then 
we may have enhanced angular correlations for massive Proca fields emitted from the black 
hole, since higher £ partial waves are less uniform. 

Another outstanding point is that for large mass, when n = 2, 3, it can be seen that the 
transmission factor starts from a constant non-zero value at the threshold oj = M {k = 0), 
at least for small £. We have checked that this does not happen for n > 4 for masses as 
large as M = 10 ~ 15, where the curves always asymptote smoothly to zero at A; = 0. Note 
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Figure 1. Transmission factors: The three rows of panels, show the first few partial waves con- 
tributing to the Hawking spectrum. Each row corresponds to a fixed mass and each column to a 
fixed dimension. In particular, the first row shows the small mass limit of the Proca theory in order 
to compare it with Maxwell's theory. 

that the parameter in the radial equations is so these are very large masses. A possible 
explanation for this phenomenon can be motivated from considerations about the range 
of the gravitational field in Rutherford scattering. In n = 2, the total cross-section for 
Rutherford scattering diverges, so the Newtonian gravitational potential is long ranged. 
This means that the effective size of the gravitational potential is infinite. The same 
happens in n = 3 but only at zero momentum k = 0. This indicates that a possible reason 
is that an incident wave at infinity with a very small momentum will still be sufficiently 
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Figure 2. Number fluxes for M = (top panels) and M — > (bottom panels): The red solid curve 
of the top panels shows the Hawking flux of particles summed over the dominant partial waves for 
the Maxwell theory. The different partial waves are multiplied by the corresponding degeneracies. 
In the bottom panels the small (but non zero) M limit of the Proca theory is shown for comparison. 
The © symbol denotes the addition of modes which are numerically equal. 

attracted by the gravitational field so that a constant non-zero fraction is still absorbed by 
the potential. In particular we note that some of the radial equations are similar in form to 
those obeyed by massive scalar and massive fermion fields, so the same effect exists for such 
fields. To our knowledge, this feature has not been noted or discussed in the literature. 
The only exception is the paper by Nakamura and Sato [35] in four dimensions, where it 
is claimed that the reflection factor for a scalar field always goes to I at a; = M (and thus 
the transmission factor goes to zero). Their result seems, however, inconsistent with Fig. 
1,2 and 3 of the paper by Page [36] (also in four dimensions), where the Hawking fluxes for 
massive fermions become constant at the k = threshold (in agreement with our result). 

Once we obtain the transmission factors, the computation of the Hawking fluxes (4.15) 
follows straightforwardly by summing up partial waves with the appropriate degeneracy 
factors (4.16). We have chosen to show the flux of number of particles. The flux of energy 
has similar features and is simply related by multiplying each point in the plots by u. 

In Fig. 2 we compare the Hawking fluxes of the Maxwell theory with the small mass 
limit of the Proca theory. For the particular case of n = 2 we have reproduced the results 
by Page [37] for the electromagnetic field and found very good agreement. All panels 
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Figure 3. Number fluxes for various M and n: (Left panel) Variation of the flux of particles for 
fixed n = 2 and variable mass. (Middle and right panels) Variation of the flux with n in a linear 
and logarithmic scale respectively. The logarithmic scale shows more clearly that the limiting flux 
at fc = is finite for n = 2,3. 

show a red solid curve corresponding to the total Hawking flux summed up over partial 
waves. The partial waves included in the sum are also represented, scaled up by the 
appropriate degeneracy factor. As claimed in the discussion of the transmission factors, as 
we increase n, partial waves with larger i become more important, for both Maxwell and 
Proca fields. One can clearly see that there is a large contribution to the total flux from the 
longitudinal degrees of freedom, since the vertical scales are larger for the Proca field. In 
particular the i = mode enhances the spectrum greatly at small energies. Note that these 
extra contributions associated with the longitudinal degrees of freedom cannot in general 
(for arbitrary mass) be described by a scalar field, since there is always a contribution 
from the coupled modes ^1,^2- That is, however, the approximation done so far in black 
hole event generators, where the W and Z fields Hawking spectra in use are those of the 
electromagnetic field (for transverse polarizations) and a scalar field (for the longitudinal 
polarization). Thus, our results can be readily applied to improve this phenomenological 
modeling. 

In Fig. 3 we show the variation of the total number flux with n and M. The left panel 
shows the expected variation with M: that the flux not only gets cutoff at the energy 
threshold u = M, but it is also suppressed with M (the same holds for n > 2). This 
is the same behaviour as found in [30, 31]. As pointed out already, in event generators 
massive vector particles are modeled using the Hawking fluxes for the Maxwell fleld and 
a massless scalar, with a cutoff at the mass threshold. In [30, 31] it was shown that 
simply imposing a sharp cut-off on the fluxes of massless scalars and fermions over-shoots 
the real amount of Hawking radiation emitted in the massive scalar and fermion channel. 
Qualitative inspection of our results suggests a similar effect for the W and Z channels in 
the evaporation. A quantitative comparison, however, requires a consideration of a Proca 
field confined to a thin brane. The middle and right panel show variation with n. In 
addition to the well known large scaling of the area under the curve and the shift of the 
spectrum to larger energies, we can also see that more partial waves start contributing 
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to the shape of the curve which becomes more wavy. This is particularly true because 
the degeneracy factors for fixed i increase rapidly with n, which is a consequence of the 
larger number of polarizations available for a vector boson in higher dimensions. Finally, 
regarding n = 2, 3 we confirm the feature that the fiux becomes a constant at = 0. This 
can be seen more clearly in the right panel in a logarithmic scale where the lines for n > 4 
curve very rapidly around that point, whereas for = 2, 3 they tend to a constant. 

6 Conclusion 

In this paper we have used a numerical strategy to solve the coupled wave equations 
obtained in the study of a Proca field in the background of a L)-dimensional Schwarzschild 
black hole. Our results show some expected features, such as the mass suppression of the 
Hawking fiuxes as the Proca mass is increased, but also some novel features, such as the 
non-zero limit of the transmission factor, for vanishing spatial momentum, in n = 2, 3. 
Moreover, for the first time a precise study of the longitudinal degrees of freedom was 
carried out. 

One application of our results will be to improve the model in the CHARYBDIS2 Monte 
Carlo event generator [22] currently in use by the ATLAS and CMS experiments. This 
simulates the production and decay of higher dimensional black holes in parton-parton 
collisions, a scenario to be soon tested at the LHC 14TeV centre of mass energy collisions. 
The Hawking evaporation can still be improved greatly through the numerical study of 
various wave equations in black hole backgrounds, which approximate the ones that could 
be produced at the LHC. This is illustrated by our results in this paper, which show the 
inaccuracy of the naive modeling of longitudinal modes by scalar fields, an approximation 
that has been used in event generators. For the phenomenology of TeV gravity scenarios, 
an important extension of the results herein is to consider a Proca field on the brane for 
scenarios where the vector bosons are confined within a thin brane. We hope to report on 
such results soon. 

Acknowledgments 

We thank Vitor Cardoso for comments on a draft of this paper. MS and MW are supported 
by FCT- Portugal through grants SFRH/BPD/69971/2010 and a research grant from the 
project CERN/FP/116341/2010, respectively. This work was partly supported by the 
project PTDC/FIS/098962/2008. 

A Functions and matrices 



2n-l 



M(r) = ^ amy' 



m 



— J. ^j.' 



.n-1 




m=0 



2n-2 



N{r) = Yl ^^y' 



m 



= n {r 



.n-1 




■m=0 



-16- 



P(r 
Q{r 
M{r 
N{r 
P{r 
Q{r 



2n~l 

m=0 
2n-2 



^a™2;"^ = zu;r"-i(2r"-i-n-l) , 

m=0 
2n 



1)' 



m=0 
2n-l 

J]/3^y'" = (n-2)r(r"-i-l)', 

m=0 
2n 

2 , j,j-2„2\ IJln-2 , , ,2„2n 



m=0 
n 

m=0 



{4 + MV2) (r2"-2 _ r"-i) + a;V2" - (n - 2) (r""! - l) 



The recurrence relations are 



- 1)03 + pI32 + 71 + (Ji) fo + o-QZ^i 



p{p + 1)02 +70 



_ a;^ + (n-l)^(p + j)(p + j-l) ia^(n - 1) ~ 



uj'^ + {n-lf{p + j){p + j-l)~ iuj{n-l) 



fi + 



^3 



(A.l) 



with 



Dj = {n-lfuj^+{uj^ + {n-l)\p + j){p + j-l)) , 

/j = - E I (am+2(p + i - m)(p + j - m - 1) + /3m+l(p + j - m) + Jrn)Pj-m + <7mZ^i-m 
m=l 

fj = - [{am+2{p + j - m){p + j - m - 1) + ^m.+l{P + j - m) + ^m)j^j-m + CFm.pj-m 



The coefficients used in the asymptotic expansion in the text are 
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The matrices used in tlie text are as follows 
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and the 2-vectors 
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